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Abstract. We investigate entanglement between electrons in serially coupled dou- 
ble quantum dots attached to non interacting leads. In addition to local repulsion 
we consider the influence of capacitive inter-dot interaction. We show how the 
competition between extended Kondo and local singlet phases determines the 
ground state and thereby the entanglement. 

1 Introduction 

In early days of quantum mechanics the entanglement between particles was considered a 
paradox. Today it has become appreciated that the ability to establish entanglement between 
qubits in a controlled manner is a crucial ingredient of any quantum information processing 
system. The interest in such systems is spurred on also by the fact that if a quantum computer 
were built, it would be capable of tasks impracticable in classical computing [T] as are, e.g., 
factoring and searching algorithms [2]. 

In general, it is desirable that the quantum computing hardware meets several criteria as 
originally proposed by DiVincenzo [3] and include (i) well defined qubits with the feasibility to 
scale up in number; (ii) the possibility to initialize and manipulate qubit states; (iii) decoherence 
processes should be minimal that quantum error correction techniques can be applied; and (iv) 
ability of detecting final qubit states as the outcome of quantum computation. It seems that 
these criteria for scalable qubits can be met in structures consisting of coupled quantum dots 
[4 1 5] which are therefore considered for implementation of quantum computing processes in 
solid state. 

In particular, recent experiments on semiconductor double quantum dot (DQD) devices have 
shown the evidence of spin entangled states in GaAs based heterostuctures [6] . It was shown 
that vertical-lateral double quantum dots may be useful for implementing two-electron spin 
entanglement [7] and it was demonstrated that coherent manipulation and projective readout 
is possible in double quantum dot systems [8]. The ability to precisely control the number of 
electrons by surface gates was also reported [§]■ 

One of the central issues regarding two-qubit operations as the basis for quantum computing 
algorithms is the creation and the control of qubit pair entanglement in a computing device pQ . 
The interaction of qubit pairs with their environment is in general a complicated many-body 
process and its understanding is crucial for experimental solid state realization of qubits in 
single and double quantum dots [5]. 

Specifically, the Kondo effect was found to play an important role in single [lOj and double 
quantum dot |llj|12|13j systems and here we report how the Kondo interaction diminishes the 
entanglement between qubits defined in DQDs even when other sources of decoherence {e.g. 
phonons) are absent. 
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Fig. 1. Schematic picture of serial DQD coupled to leads. 



2 Coupled quantum dots with interaction 

We consider a serially coupled DQD: a device with the ability to produce entangled pairs that 
may be extracted using a single-electron turnstile [13] ■ We model DQD using the two-impurity 
Anderson Hamiltonian 



where c is creates an electron with spin s in the dot i — A or i — B and rii S = c is Ci S is the 
number operator. The on-site energies e and the Hubbard repulsion U are taken equal for both 
dots. The dots are coupled to the left and right noninteracting tight-binding leads with the 
chemical potential set to the middle of the band of width 4fo. Each of the dots is coupled to 
the adjacent lead by hopping t' and the corresponding hybridization width is f = (i') 2 /io- 
Schematically this setup is presented in Fig. [I] The dots are coupled capacitively by a inter-dot 
repulsion term Vn A n B . 

In this paper we concentrate on the low temperature properties of DQD system determined 
from the ground state. We expand the ground state in the Schonhammer and Gunnarsson 
projection-operator basis |15|16j |#av) = P{\i)P(\>j) |0), which consists of projectors P\i where 
i £ {A,B}- e.g., P (0l) = (1 - n^) (1 - n^), P {ll) = J2a n ™ i 1 ~ n «<?)> p (?i) = n tf n H ~ and 
additional operators involving the operators in leads. We used up to ~ 100 additional com- 
binations of operators consisting of, for example, Pi^zji) — P(oi)VjP(u)j where vj denotes the 
tunneling to/from dot i to the site j in the lead. These operators are applied to the state |5), 
which is the ground state of the auxiliary noninteracting DQD Hamiltonian of the same form 
as H 7 but with U, V = Q, renormalized parameters e, t, t' — > e, t, t' and additional parameter t" 
which corresponds to hopping from left dot to right lead and vice versa which although absent 
in the original Hamiltonian is present in the effective Hamiltonian in some parameter regimes. 

The starting point towards the understanding of the ground state of DQDs are the filling 
properties of isolated DQDs. The first electron is added when e = t, and the second when 
e = -t + J - {{U + V) - \U - V\}/2, where J = [-\U - V\ + y/(U - V) 2 + 16i 2 ]/2 is the 
difference between singlet and triplet energies. For e + U/2 + V = Q DQD is doubly occupied, 
n = (n A + n B ) = 2, and the ground state is ^=[a(|tl) - lit)) + (3(\20) - |02})|, where a//3 = 

it/(V — U+ y/(U — V) 2 + 16i 2 ). Here we use notation |ti) corresponding to spin-up and spin- 
down states on sites A and B, and 1 20) to double and zero occupancy of sites A, B. The range 
of e where single occupation is favorable is progressively diminished when V ^ U. For large t 
or at (and near) V — U the molecular bonding and anti-bonding orbitals are formed as is seen 
here from a ~ (3. 

When DQDs are attached to the leads the low temperature physics is to the large extent 
the same as that of the two-impurity Kondo problem studied by Jones, Varma and Wilkins 
two decades ago |17|18j . There two impurities form either two Kondo singlets with delocalized 
electrons or bind into a local spin-singlet state which is virtually decoupled from delocalized 
electrons. The crossover between the regimes is determined by the relative values of the ex- 
change magnetic energy J and twice the Kondo condensation energy, of order the Kondo tem- 
perature given by the Haldane formula, Tk = y/U r/2 exp(— 7re(e + U)/2r). Such results were 
obtained by the analysis of a two-impurity Anderson model by means of slave-boson formal- 
ism |19I20I21I22I23| . numerical renormalization group (24125126] or present formalism [26127] . 



H= ( en i + Un^ni\) + Vn A n B - t ^(c 
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Resembling behavior was found also in particular regimes of triple quantum dot systems .28], 
and DQDs in side coupled [29] and parallel [30] configurations. 

3 Entanglement 
3.1 Spin entanglement 

Quantum entanglement as a physical resource was first defined for two distinguishable particles 
in a pure state through von Neuman entropy and concurrence [31.32 33.34J. However, amongst 
the realistic systems of major physical interest, electron-qubits have the potential for a much 
richer variety of entanglement measure choices due to both their charge and spin degrees of 
freedom. In systems of identical particles, for example, generalizations are needed to define 
an appropriate entanglement measure which adequately deals with multiple occupancy states 
[35I36I37I38| . 

When entanglement is quantified in fermionic systems the measure must also account for 
the effect of exchange [39] as well as of mutual electron repulsion. In lattice fermion models 
entanglement is sensitive to the interplay between charge hopping and the avoidance of double 
occupancy due to the Hubbard repulsion, which results in an effective Heisenberg interac- 
tion between adjacent spins [3(1]. Entangled fermionic qubits can be created with electron-hole 
pairs in a Fermi sea [41] and in the scattering of two distinguishable particles [42] . A spin- 
independent scheme for detecting orbital entanglement of two-quasiparticle excitations of a 
mesoscopic normal-superconductor system was also proposed recently [43j . 

For two distinguishable particles A and B, described with single spin-i (or pseudo spin) 
states s =| or j and in a pure state <Pab) = J2 SS ' ass ' I s )a|s')b concurrence as a measure of 
entanglement is given by [32] 

Co = 2|a TT a ii - a u a a \. (2) 

Two qubits are completely entangled, C = 1, if they are in one of the Bell states [31] , e.g., 
singlet |^ab) cx | tl) - | IT). 

A qubit pair represented by two electrons in DQDs and in the contact with the leads 
acting as a fermionic bath can not be described by a pure state and entanglement can not 
be related to the concurrence given with the Wootters formula Eq. ([2]). In the case of mixed 
states describing qubit pairs concurrence is related to the reduced density matrix of the DQD 
subsystem 34 44 45 , where for systems that are axially symmetric in spin space the concurrence 
may conveniently be given in the closed form [46j , 

C = max(0, C^, Cm), 

C n = 2\(S+S B }\ - 2yJ{PlPl){PiP&, (3) 
C ll =2\(S+S+)\-2y/(PlPi)(PiPl), 

where = (S^Y = c^Cii is the electron spin raising operator for dot i — A or B and 
Pf = rii s (l — ni,-s) is the projection operator onto the subspace where dot i is occupied by one 
electron with spin s. 

In the derivation of concurrence formula Eq. ([3]) the reduced density matrix was obtained 
by projecting onto four local spin states of | f )a, | I )a, | T )b, and | J. )b, corresponding to singly 
occupied DQD sites A and B, respectively. If t/U is not small the electrons tunnel between 
the dots and charge fluctuations introduce additional states with zero or double occupancy of 
individual dots [35140] . As pointed out by Zanardi 40 in the case of simple Hubbard dimer the 
entanglement is not related only to spin but also to charge degrees of freedom which emerge 
when repulsion between electrons is weak or moderate. 

For systems with strong electron-electron repulsion, charge fluctuations are suppressed and 
the states with single occupancy - the spin-qubits - dominate: the concept of spin-entanglement 
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quantified with concurrence can still be applied. We use spin-projected density matrix and 
consider only entanglement corresponding to spin degrees of freedom. Due to doubly (or zero) 
occupied states arising from charge fluctuation on the dots (caused by tunneling between the 
dots A and B or due to the exchange with the electrons in the leads), the reduced density 
matrix has to be renormalized. The probability that at the measurement of entanglement there 
is precisely one electron on each of the dots is less than unity, Pn < 1, and the spin-concurrence 
is then given with 

C = C /P n> (4) 

where P n = P n + P,| , and P n = (P^ + P A P^), p \\ = ( p a p b + p a p b) are probabilities for 
antiparallcl and parallel spin alignment, respectively. Such procedure corresponds to the mea- 
surement apparatus which would only discern spins and ignore all cases whenever no electron, 
or a electron pair would appear at one of the detectors at sites A or B. 

3.2 Charge (isospin) entanglement 

In the ground state of two isolated impurities coupled by a capacitive (but not tunneling) term 
V = U, 4 'spin states' \a1a2) and 2 'charge states' |20), 1 02) are degenerate. By introducing the 
pseudospin operator [47] T % = 1/2 Ylw=i 2 J2a c i<r T u ,c l'<n where t 1 are the Pauli matrices, and 
the combined spin-pseudospin operators = S i T^ , the Hamiltonian is SU(4) symmetric. 
As long as the SU(4) symmetry breaking terms are small enough {e.g., tunneling t — > 0) the 
ground state of such DQDs attached to the leads remains close to an SU(4) state with 'spin' 
screened by the electrons in the leads [27] . 

If V » U charge states dominate and in this case charge concurrence can be defined in a 
direct analogy with the previous spin case. In Eq. ([3|) one just has to replace the spin operators 
with their corresponding isospin counterparts, e.g., = c^c^i = (S^Y — > = c\^c\i = 
(T+)t for sites A =A,B and S*| = (n AT - n Xi )/2 -> T{ = (n x - l)/2. If the probability for spin 
states is significant, appropriate renormalization to charge states is analogous to Eq. (j4|), but 
with corresponding isospin operators. The density matrix is here renormalized with the proba- 
bility P20 that precisely ua,b = 2 and ne,A = electrons occupy individual dots (corresponding 
to apparatus which only measures occurrence or absence of pairs at each of its detectors) . 



4 Numerical results 

4.1 Concurrence 

We here present the results for zero temperature (ground state) concurrence of a qubit pair in 
DQDs in the absence of magnetic field. Temperature dependence of concurrence for the case 
of V = is given in Ref. [25] . Expectation values (...) in the concurrence formula Eq. ((3]) are 
now calculated using the ground state therefore {S^S^} = and Cm < 0. We consider the 
particle-hole symmetry point with n — 2 and e + U/2 + V = 0. 

Qualitatively, the concurrence is significant whenever enhanced spin-spin correlations in- 
dicate inter-dot singlet formation. As shown in Fig. [2](a) for UjP = 12 and r/to = 0.1, the 
correlation function (Sa ■ Sb) tends to —3/4 for J large enough to suppress the formation of 
Kondo singlets, but still J/U < 1, that local charge fluctuations are sufficiently suppressed. In 
particular, the local dot-dot singlet is formed whenever singlet-triplet splitting superexchangc 
energy J > J c ~ 2Tk- With increasing V — ► U, and above U, the probability for singly occupied 
spin states, Pn — Pf| + Pi is significantly reduced, Fig. H£b) and Fig. [3] which also leads to 
reduced spin-spin corelation, Fig.[2ja). In this limit the concept involving isospin entanglement 
can be applied (not shown here). 

Concurrence, corresponding to the correlation function from Figs. [2ja,b) is presented in 
Fig.[3]for various values of V. As discussed above, C is zero for J below ~ 2Tk due to the Kondo 
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Fig. 2. (a) Spin-spin correlation for V/U = 0,1/3,2/3,1 (full lines) and V/U = 5/4 (dotted) for 
U/r — 12,r/to = 0.1. (b) Probabilities for parallel (lower curves - dashed) and anti-parallel (upper 
curves - full) spins of electrons in the DQD for V/U ratios as in (a). Note that the probability for parallel 
spins for V/U = 5/4 is almost zero (dashed-dotted) , while P n < 1/2 for J/T K < 1000 (dotted); the 
probabilities do not sum to 1. The deficiency (which goes to zero as U — > oo) is due to states with 
double particle (or hole) occupancy on at least one dot. 




j ' i 



1000 



Fig. 3. Concurrence (full curves) and single particular dot occupation probability Pn (dashed) for 
V/U = 0, 1/3, 2/3, 1. For V/U = 5/4 the concurrence and Pn are plotted dotted and dashed-dotted, 
respectively. Parameters are as in Fig. [2] 



effect, which leads to entanglement between localized and conducting electrons [48J instead of 
the A-B qubit pair entanglement. In finite magnetic field irrespectively of temperature the 
concurrence abruptly tends to zero for B > J (not shown here) [49) . 



4.2 Conductance 



One of the most directly measurable properties of DQDs is the linear conductance. We calculate 
the zero temperature conductance using the sine formula (SF) [50151152) . G = Gosin 2 [(£ + - 
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Fig. 4. Conductance of DQD as a function of gate voltage and inter-dot tunneling rate for U/T = 
20,r/to — 0.04. Pictograms indicate dominant ground state regimes: molecular orbital Kondo effect, 
local spin-singlet formation and two separate Kondo effects. 

E_)/it()L], where Go = 2e 2 /h and E± are the ground state energies of a large auxiliary ring 
consisting of L non-interacting sites and an embedded DQD, with periodic and anti-periodic 
boundary conditions, respectively. Alternatively conductance can be obtained also from the 
Green's function (GF) corresponding to the effective noninteracting Hamiltonian H with ef- 
fective parameters [50153) . The advantage of the former method is better convergence in the 
strong coupling regime, however, the accuracy of the SF method depends only on the accuracy 
of the ground-state energy and is therefore in some cases more robust. By comparing results of 
both methods we checked for the consistency and the convergence. 

The conductance as a function of interdot hoping t and e + U/2 in the absence of inter- 
dot repulsion, V = 0, is presented in Fig. [4] The Hubbard repulsion is set to U/T — 20 and 
hibridization to r/to = 0.04. As in the previous Section three different regimes of i correspond- 
ingly reflect in the results for conductance. For large t/T > 1 (but with r/U « 1) the DQD 
is in molecular-orbital Kondo regime when occupancy is odd, i.e, n ~ 1, 3. Typical Kondo con- 
ductance plateau of width ~ U/2 is developed around bonding (and anti-bonding) molecular 
orbital level it. These two unitary conductance regions with reducing t become progressively 
sharper when we enter t/U <K 1 regime. There the description of DQD in terms of bonding/anti- 
bonding orbitals should be replaced with local picture. Due to strong electron-electron repulsion 
local charge fluctuations are suppressed at the point of particle-hole symmetry with n = 2 and 
e + U/2 = 0. Thick full line corresponds to points of G = Go and there exists some critical t c 
where two conductance peaks merge (bullet) and for t <t c conductance is less then Go (dashed 
line). The corresponding critical superexchange interaction is of the order of Kondo tempera- 
ture as before, J c ~ 2Tk- As discussed above, in this regime each of the dots undergoes the 
Kondo effect where local moment is screened by conducting electrons in the adjacent lead and 
left - right sides of the system become decoupled which leads to vanishing A^B conductance 
as G oc {t/rf [19]. 

5 Summary 

The main results concerning entanglement of qubit pairs in serially coupled double quantum 
dots are extracted in Fig. The charge fluctuations An A = (nj.) — (ua) 2 , contour plot in 
Figs. [5ja), are suppressed for sufficiently large repulsion, e.g., U/T > 10. In this limit and in 
vanishing magnetic field, the DQD can be described in terms of the Werner states [54] and 
becomes similar to recently studied problem of entanglement of two Kondo spin impurities 
embedded in a conduction band [55]. In this case, C^i ~ 2(— (Sa • Sb) — j) ~ P\i — 2P\\ for 
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U/T V/U 

Fig. 5. (a) Charge fluctuations (short-dashed), Sa • Sb = —1/4 (long-dashed) and C = (full) in the 
(U/T, J/T K ) plane, (b) S A ■ S B = -1/4 (long-dashed) and C = (full) in the (V/E7, J/Trt) plane. P n 
is shown with dotted line (right scale). 

C|| > 0. For large C//-T, where the charge fluctuations vanish, the (Sa • Sb) = — h line (dashed 
line) progressively merges with the C = boundary line (full) . 

In Fig. H^b) phase diagram with fixed U/T — 12 and corresponding to Fig. [5(a) presents 
V/U dependence of C = boundary line (full) in comparison with the (Sa ■ Sb) = —\ line 
(dashed). With V exceeding U the probability for well defined spin-qubit pairs in DQD, Pn, 
rapidly decreases which means that states with doubly occupied or empty individual dots dom- 
inate. In this regime C = line is pushed to much lower J/Tk- For V > U the probability Pn 
becomes progressively negligible giving more meaning to considering charge (isospin) entangle- 
ment instead. It should be noted, however, that in realistic DQD systems intersite repulsion V 
is in general weaker compared to U and that this regime would not be easily reached exper- 
imentally One possibility, where ^-interaction could dominate, are systems with strong local 
electron-phonon interactions which may significantly renormalize local U [56] without affecting 
capacitive interaction V . 

To conclude, we have found generic behavior of spin-entanglement of an electron pair in 
serially coupled double quantum dots. On the one hand, we have shown quantitatively that 
making the spin-spin exchange coupling J large by increasing tunneling t, leads to enhanced 
charge fluctuations, whilst on the other, at small magnetic interactions J < J c entanglement is 
suppressed as the DQD system undergoes the Kondo effect. Various regimes are explained and 
supported with typical numerical examples. 

We thank T. Rejec for his GS code and suggestions. We acknowledge J. Bonca and R. Zitko 
for useful discussions. We also acknowledge support from the Slovenian Research Agency under 
contract Pl-0044. 
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